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I. INTRODUCTION 

Electroweak processes are an important tool in under- 
standing the Standard Model (SM) of particle physics, 
serving as an input into tests of the unitarity of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix and as a 
probe for new physics. The hadronic matrix elements 
that characterise the strong interaction dynamics of these 
processes are a crucial ingredient in the determination of 
CKM unitarity. 

Global fits to the CKM unitarity have, in recent years, 
indicated some tensions at the 2-3(7 level within the SM 
[T}]3]. In many cases, the constraints on the CKM uni- 
tarity triangle are limited by the precision with which 
the nonperturbative inputs are known and thus it is im- 
perative that these inputs are determined as precisely as 
possible. 

The HPQCD collaboration has undertaken a suite of 
precision calculations of heavy-light mesons as part of 
a program to precisely determine nonperturbative con- 
tributions to electroweak parameters. Recent calcula- 
tions of the decay constants Jb and fs, have achieved 
a precision at the 2% level, by taking advantage of the 
small discretisation errors and good chiral properties of 
the Highly Improved Staggered Quark (HISQ) action 
[3 [6] . These results represent the most precise currently 
available for these decay constants. In addition, nonper- 
turbative studies of the heavy-light semileptonic decays 
B TTliy, B K£+£~ and Bg Kiiy are underway [YJ. 

The work of Ref. [5] and [7j use HISQ light quarks 
and the nonrelativistic QCD (NRQCD) action for the 
heavy quarks. These calculations require matching the 
heavy-light axial-vector and vector currents in the effec- 
tive theory on the lattice with full QCD. In this article 
we report on the one loop perturbativc matching of the 
HISQ-NRQCD axial-vector and vector current matching 
for both massless and massive HISQ quarks. As part of 



this procedure we determine the mass and wavefunction 
renormalisation for massive HISQ quarks. Our match- 
ing results for massive HISQ quarks will be relevant for 
future studies of heavy-heavy decays — >■ 

In the next section we describe the quark and gluon 
actions used in our calculation. We then review the for- 
malism for extracting renormalisation parameters from 
relativistic lattice actions and apply these procedures to 
first massless and then massive HISQ quarks. We include 
results for the one loop NRQCD mass and wavefunction 
renormalisation in Section llll El In Section HVl we outline 
the calculation of the matching coefficients and then, in 
Section |Vj we present our results for a range of heavy 
quark masses. We conclude with a summary in Section 

ED 



II. THE LATTICE ACTIONS 

A. Gluon Action 

We use the Symanzik improved gluon action with tree 
level coefficients [8l411j. given by 



So = 



3lin 



E 

X.fJ,>V 



(1) 



Here P^i, is the plaquette, 

P^.. = ^RcTv{u^{x)U,{x + fj)Ul{x + J>)t/^(.T)}, (2) 
and R^u the six-link loop, 

Ruu = ^ReTr\u^{x)U^,{x + fi)U,{x + 2fi) 

Ulix + fi + mli^ + i^Mix)], (3) 



X 
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TABLE I. Summary of quark mass notation. 



HISQ 


amo 
amtrcc 

ami 
amkin 


bare light quark mass 
tree level pole mass 
one loop pole mass 
kinetic mass 


NRQCD 


aMo 


bare heavy quark mass 



with f3 = 2Nc/g'^ and uq the tadpole improvement fac- 
tor [12j . Radiative improvements to the gluon action do 
not contribute to the one loop matching calculation. In 
general, radiative improvement generates an 0{as) in- 
sertion in the gluon propagator. There are no external 
gluons in our calculation, so any such improvements only 
contribute at two loops and higher. 
We include a gauge-fixing term 



H 2 



(4) 



where is the symmetrised difference operator, which 
acts on the gauge fields as 



A„ 



A,, 



(5) 



and ^ is the gauge parameter. Where possible, we confirm 
that gauge invariant quantities are independent of the 
choice of gauge parameter by working in both Feynman, 
^ = 1, and Landau, ^ = 0, gauges. 



B. Light Quark Action 

We discretise the light quarks in this work using the 
Highly Improved Staggered Quark (HISQ) action ^15] . 
The HISQ action significantly reduces taste breaking dis- 
cretization errors and has been used successfully to sim- 
ulate both b and c quark systems [TH - fTT] . There are 
two equivalent methods for writing staggered quark ac- 
tions, using either four component "naive" fermions or 
one component "staggered" fields [TSl [TU] . Throughout 
this calculation we use the naive fermion representation 
and we denote the bare quark mass amg. In Section 
|IIID 1 1 we present our results for massless HISQ quarks, 
corresponding to aniQ = 0. Before we present the quark 
actions used in this work, we pause to briefly discuss 
some notation, which we summarise in Table |lj We use 
four different quark mass definitions for relativistic HISQ 
quarks: the bare quark mass; the tree level and one loop 
pole masses, amtree and ami respectively; and the ki- 
netic mass, amkin- We distinguish these relativistic quark 
masses from the nonrelativistic quark mass in NRQCD 
by using a lowercase m for HISQ quarks and an upper- 
case M for NRQCD quarks. Only the bare heavy quark 



mass aMo is required for nonrelativistic quarks in this 
calculation. 

The starting point for constructing the HISQ action is 
the AsqTad action [18 , which is given by 

^AsqTad =a^J2 ^(^) ht^^t'""^"" + ™0) ^X), (6) 

X 

where the AsqTad operator is 



7AsqTad _ - — CV \^ 
^l ~ ^ 6 ^ 



(7) 



Here the three-link term (V^)'^ is referred to as the 
"Naik" term and the superscript F indicates that we use 
fattened links in the lattice difference operator V^. The 
fattened links are given by 



where 



X- AsqTad 




(8) 



(9) 



(10) 



symmetrised 



The second term in Equation ([9]) is the so-called "Lep- 
age" term. The difference operator acts on fermion fields 
as 

^M^) - ^ {u^,{x)^^{x + fi)- ulix - A)^(x - A)] , 
(11) 

whilst the discretised derivatives acting on link variables 
are, for n ^ v, 



V^C/^(a;) 



V(^)C/.(2;) 
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U^ix)U,ix + fi)Ulix + 0) 

ulix - fi)U^{x - p.)Uf_,{x - ft + v) 



(12) 



U^ix)U,ix + (i)Ul{x + 2U,ix) 
+ U^ix - jji)Ui,{x - jjL)U^{x - jl + v) 



(13) 

The HISQ action is an extension of the AsqTad action 
that includes two levels of link fattening and a tuned co- 
efficient for the "Naik" term. Whilst the AsqTad action 
has negligible tree level errors for light quarks, this is not 
true for charm or bottom quarks [f3ji. Charm quarks are 
generally nonrelativistic in typical mesons, so the rest 
energy of the quark is much larger than its momentum. 
The dominant tree level errors are therefore ©(a^mp). 
One suppresses these errors by tuning the coefficient of 
the Naik term 



a 
"6 



a 



(l + f)(V^)^ 



(14) 
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One also adds a second level of fattening in the link 
variables to reduce the discretisation errors arising from 
taste exchange interations in the HISQ action. Between 
the smearing operations, one sandwiches a reunitarisa- 
tion operator, lA^ that projects the smeared link variables 
back to SU{2>) or C/(3). For simplicity, the Lepage term 
is included in the HISQ action only after the second level 
of link fattening. The resulting action is 

5hisq = a'E^(^) (7pV™SQ ^.(a;), (15) 



and the correction terms are 



where 



{UF) 



(16) 



The superscripts indicate that the first operator, 
is built from the full HISQ-smeared links, given by 



«'(Vp)2 



(17) 



whilst the second operator, V^^^\ uses only one level of 
smearing: 



j{UF) _ 



(18) 



We define the operator in Equation (10). 



We give results for both massless and massive HISQ 
quarks. For massless quarks the tuning parameter is just 
e = 0. For massive quarks we set the tuning parameter 
to its tree level value, e = etrco, for consistency with non- 
perturbative simulations [6j. We discuss this in more 
detail in Section UlIDI 



C. Heavy Quark Action 

For the heavy quark fields, V'lx, i), we use the NRQCD 
action of PUI [3T], which is improved through 0(1/Mg) 
and 0{a^) and includes the leading relativistic (^(I/Mq) 
correction. The NRQCD action is 



'S'nrqcd = E -0* V't - "0* ( 1 - I 1 - 



yi,t 



X (7] ( 1 



2n 



2 J \ 2n 
aSH 



1 



^t-i, (19) 



where ipl — ip^(x.,t) and ^t_i = 'ijj{x.,t — 1). 

Here the leading kinetic term in the NRQCD action is 
given by 



A(2) 
' 2aMn ' 



(20) 



aSH = — ci 



C3 



C4 



(A( 



C2: 



8(aA//o)3 8(aMo)2 
1 



V-E-E- V 



8(aMo)2 
1 



cr- VxE-ExV 



2aMr 



CT • B + C5 — — — C6 



24aMo 16n(aMo)2 ' 
(21) 



All the derivatives are tadpole improved and the discre- 
tised difference operators are 



V 



(3) 



(22) 

where the improved operators act on fcrmion fields via 

Vf V(a:) = C/,,(a;)V^(.T + A) 

+ Ul{x-fi)i^{x-fi)-2i;{x), (23) 
1 ' 



Ul{x~fi)Ul{x-2^)^:{x-2fi) 
Uf,{x)ip{x + fi) + Ul{x- fi)il;{x - jj.), 



(24) 

V^^^ix) = U^{x)U^{x + fi)^{x + 2A), 

+ Ul(x- [i)Ul{x - 2fi)ij{x ~ 2A) + 6i^{x) 

- 4(Uf,ix)il;{x + fi) + Ulix - fL)^ix - A)) . 

(25) 

The improved chromo-electric and -magnetic fields, Ej = 
and Bj = —eijkFjk/2, are defined in terms of the 
improved field strength tensor, given by [19] : 

- If^A^) - I (u^{x)F^,{x + fi)Ul{x) 

+ Ul{x - il)Ff,^{x - p.)Uf,{x - A) - (m ^ i^)^ , 

(26) 

where 

F^,ix)= -^^{n,Ux)-nl^ix)), (27) 
fi^^(x) = ^ Ua{x)Ui3{x + a) 

{(",/3)} 

xU^^{x + a + p)U-p{x + p). (28) 
The final sum runs over 

{{a, 13)} = {(/i, i^), (i^, -m), (-M, ->^), {-1^, m)}, (29) 
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with y^^v. 

The values of the coefficients, c^, in the NRQCD action 
are fixed by matching lattice NRQCD to full QCD. We 
use the tree level values oi ci — \ for all i = 1, • • • , 6, 
and do not consider the effects of radiative improvement 
of the NRQCD action. 



III. QUARK SELF ENERGY 

Perturbative calculations of the self energy for massless 
AsqTad quarks were carried out in [22] as part of the 
matching calculation for NRQCD-AsqTad currents. In 
this work, we extend these results to HISQ fermions. We 
update the results for the massless case and generalise 
the results to massive quarks, applying the methods of 
[25] to extract the self energy parameters. 



A. HISQ Parameters 

The general formalism for self energy calculations is 
laid out in and developed in [23] ■ In this section we 
apply this formalism to the HISQ action, concentrating 
on the massive case. 

We start with the quark two-point correlation function. 



In Sections jlllD Ij and jlll D 2 we will use Zq and Zq to 



denote the massless and massive wavefunction renormal- 
isations respectively; in this section, however, we use Zq 
as shorthand for either Zq or Zq for notational simplicity. 

We renormalise at the point (poiP) — {iEjO) and 
therefore consider a zero spatial momentum quark prop- 
agating forward in time, for which one expects 



G{t,0) = ZQ 



-Et 



1 



- 70 



(34) 



We denote the momentum space quark propagators for 
the full and free theories G{p) and Go{p) respectively. 
These propagators are related via the quark self energy. 



(35) 



(V'(i,p')V'(0,p)) = (27r)35(p-p')G(t,p), (30) with Equation jSf I 



where the self energy is the sum of all one-particle irre- 
ducible graphs; in perturbation theory one assumes that 
the self energy is a "small" correction. The pole corre- 
sponding to the single particle quark state has a non-zero 
residue in the limit that the self energy vanishes, whilst 
the residues of the multiparticle states vanish in the ab- 
sence of an interaction. 

Carrying out the Fourier transform in po of the full 
quark propagator, G{p), one finds 

G{t,p) ^ ^e-^P"'Gipo,p). (36) 

We identify this expression at zero spatial momentum 
which enables us to relate the mass 



which defines the quark propagator G'(i, p). The bare 
quark field "0(0, p) creates multiparticle states in addi- 
tion to a one-quark state and so one expects the quark 
propagator to take the form 



G(t,p) = Z2(p)e-^(P)*rp, 



proj 



(31) 



and wavefunction renormalisation to parameters in the 
action, via the quark propagator. In the following deriva- 
tions, we will neglect factors of the lattice spacing a for 
simplicity. These can be easily included at the end of the 
derivations by dimensional analysis. 



B. Pole mass 



Here Tpi-oj is a projection operator in Dirac space; the el- 
lipses represent multiparticle states and lattice artifacts, 
which we will not consider any further; and (p) is the 
single quark residue. 

The use of a lattice regulator distorts the mass shell of 
the quark, which would otherwise satisfy the relativistic 
dispersion relation in Euclidean space. To account for 
the distorted pole position in a systematic manner, one 
therefore defines the rest mass of the quark, toq, as 



vfiQ = S(p = 0) 



and the wavefunction renormalisation as 



Zq = Z2(p = 0). 



(32) 



(33) 



For HISQ fermions, the form of the free propagator is 
^Q^^'P) ^ X! sin {p^^) Kf,{p) + viq. (37) 



Here toq is the bare quark mass and 

1 + e 2 
K^,{p) = 1 + (sinp^) , 



(38) 



We write the one loop self energy as 

Ed) (p) = ^ sin (p^) Ef (p) + eW (p)I, (39) 



where I is the identity element of the Clifford algebra, so 
that the one loop propagator is 



5 



(sin {pp)f K,{p) - asY.f{p) " + [mo - a^^ip)] 



(40) 



I 

At zero spatial momentum the pole condition for the where Zq^ is the tree level wavefunction renormalisation, 



forward propagating quark is 
1 + e 



sinh(i?) (l - ^ (sinh(i?))2 - asY^^i^^ = mo-a,Y^^\ 

(41) 

where we have neglected the arguments of E^'*^ and S'^' 
for clarity. We now expand the quark energy and tuning 
parameter e to one loop as 



given by 

^ = <^ COsh(TOtrGo) 



Z. 



1 Z (smh(mtrcc)) 



(47) 



E = rritroc + a^TOi, 

e = etreo + ttsCl- 



(42) 
(43) 



C. Wavefunction renormalisation 



We now extract the wavefunction renormalisation from 
Substituting these expressions into the pole condition, the quark propagator. Recall that the wavefunction 



Equation gives an expression for the tree level pole renormalisation is given by residue of the single particle 
mass, mtreo, at fixed bare mass, mo: momentum pole obtained by identifying Equations p4| 

1 + etroo 



sinh(TOtroc) 



1 - 



6 



(sinh(mtrco))' 



and ( 36 ) at zero spatial momentum, whence 



mo. (44) 



We then fix etroo by requiring that the tree level pole mass 
is equal to the tree level kinetic mass. We discuss this 
condition in more detail in Appendix [K\ One ultimately 
finds 



■K j a 



IT / a 



27^ 



ipot 



G{po,Q)^ZQe 



1 + 70 



(48) 



It is convenient to re-express this relation in terms of 
the variable z = 6*^" : 



1 



1 



(sinh(mtroc)) 



12mtrc 



cosh(mtroo) sinh(mtrco) 



— -K I a 



27r 



e-*P°*G(po,o; 



2tt 



(45) 



(49) 



where the contour of integration is now around the unit 
circle in the complex z-plane. Writing the propagator as 
Expanding this equation gives Equation (24) of [H]. We G(z, 0) = gi{z)/g2iz), the residue at z = zi = e"-^ is 



obtain a precise numerical value for the tree level mass 
by solving Equations ( 44 1 and ( 45 1 self-consistently; we 



find that a series solution is insufficiently accurate for our 
accurately setting the light quarks onshell. 
We repeat the process at one loop to obtain 

mi = Z^") {| (sinh(mtrec))' + sinh(mtrcc)s[,^) - S®} , 



Res {z*-iG(z,0)} 



giizi) 



(50) 



"2152(^1) 

where the prime indicates differentiation with respect to 
z. 

In this case the quark propagator is given by Equation 



(46) (40) and we obtain 



Res {z*^iG(z,0)} =e-^*i4r^J cosh(£;) 



1 - i^(sinh(£;))2 



+ asi 



d r 

dpo 



isin(po)4'''+sW 



(51) 



Comparing this equation with Equation ( 48 1 , we read off the wavefunction renormalisation as 

/■ \ -1 

d 



cosh(i?) 



1 - ^-^{smh{E)f 



zsin(po)S^^^+S™ 



(52) 



We again expand the mass and tuning parameter as in Equations ( 42 1 and ( 43 1 . The tree level result reduces to 
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Equation ( 47 1 , whilst the one loop wavefunction renormalisation is 



z9.' = Z. 



(0) J ei 
Q 



cosh(r7itrcc) (sinh(mtrcc)) - rui sinh(TOtrco) 



sin(po)S^'*^-*S« 



1 - 



1 + etrc 



2 (cosh(TOtrcc))^ + (sinh(mtrcc))^ 



(53) 



Here we have found it convenient to factor out the tree level wavefunction renormalisation to ensure the one loop 
terms have the correct infrared divergence In other words, we set 



Q 

.(0) 



l + a,Z«+0(a2). 



(54) 




FIG. 1. Contributions to the one loop self energy required 
for the HISQ quark mass and wavefunction renormalisation. 
On the left is the "rainbow" diagram and on the right the 
"tadpole" diagram. Straight lines represent light quarks and 
curly lines indicate gluons. 



D. Numerical Results 

In this section we summarise our results for both mass- 
less and massive HISQ quarks. The diagrams that con- 
tribute to the self energy at one loop are shown in Fig- 
ure ^ The continuum- like contribution is the "rainbow" 
diagram, shown on the left of Figure [l] On the right is 
the lattice artifact "tadpole" diagram. We calculated the 
corresponding Feynman integrals using two independent 
methods. 

Our first method employed the automated lattice per- 
turbation theory routines HiPPy and HPsrc [55J[5^. These 
routines have now been used in a number of perturbative 
calculations, for example in pT l [2T | [25 l [27H32] . and have 
been extensively tested against results published in the 
literature. 

Evaluating the relevant Feynman integrals is a two 
stage process: we first generate the Feynman rules with 
HiPPy, a set of Python routines that encode the Feynman 
rules in "vertex files" . These vertex files are then read in 
by the HPsRC code, which is a collection of Fortran mod- 
ules that reconstruct the diagrams and evaluates the cor- 
responding integrals numerically, using the Vegas algo- 
rithm All derivatives of the self energy are imple- 
mented analytically using the derived t ay lor type, de- 
fined as part of the Fortran TaylUR package [M]. We 
performed our calculations on the Darwin cluster at 
the Cambridge High Performance Computing Service, as 
part of the DiRAC facility, and the Sporades cluster at 
the College of William and Mary with routines adapted 
for parallel computers using MPI (Message Passing In- 



terface) . 

In contrast to previous matching calculations, such as 
[22] , we do not attempt to present Feynman rules for the 
improved NRQCD and massive HISQ actions: the au- 
tomated lattice perturbation theory procedure does not 
require such explicit analytic expressions. This method 
therefore reduces the possibility of algebraic errors in the 
manipulation of Feynman integrands. 

We undertook a number of tests of our automated per- 
turbation theory code. In particular, we reproduced the 
results of [H] with massless AsqTad light quarks and 
NRQCD heavy quarks. The chief advantage of the auto- 
mated lattice perturbation theory routines is the relative 
ease with which different actions can be implemented in 
the calculation. Once the correct HPsrc code is in place 
to calculate the requisite Feynman diagrams, switching 
actions is just a matter of replacing the input vertex files 
generated by HiPPy. 

In many cases, we established that gauge invariant 
quantities, such as the mass renormalisation, are gauge 
parameter independent by working in both Feynman and 
Landau gauges. 

Furthermore we confirmed that infrared divergent 
quantities, such as the wavefunction renormalisation, ex- 
hibited the correct continuum-like behaviour. We reg- 
ulate the infrared behaviour with a gluon mass for 24 
different values of the gluon mass between a?\^ = 10~^ 
and a^A^ = 10^^^. Fitting these results to a logarith- 
mic function establishes that the code correctly repro- 
duced the expected logarithmic behaviour. To extract 
the infrared finite piece of infrared divergent quantities 
we constrain the fit function to have the correct logarith- 
mic divergence. 

At finite lattice spacing offshell contributions to the 
vertex renormalisation must be removed to restore the 
correct continuuum-like infrared behaviour. We set the 
HISQ quarks exactly onshell and remove offshell con- 
tributions to the vertex renormalisation with an onshell 
projector. This corresponds to imposing the equation of 
motion on the quark or antiquark spinor, just as would 
be done analytically |35j . It is important to ensure the 
quark is set exactly onshell, by solving the full inverse 
tree level HISQ propagator for the timelike component 
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of the quark momentum, or the continuum infrared be- 
haviour is not recovered. We found that this requires very 
precise values for mtioc and etioo (see Table III ) and that 
using only a few digits is insufficient. Likewise the equa- 
tion of motion for the massive HISQ propagator must be 
exactly satisfied for the ofFshell contributions to be fully 
removed. 

Our second method is based on Mathematica and 
Fortran routines developed previously for matching of 
NRQCD/AsqTad currents |22] and adapted here for 
HISQ light quarks. Although Feynman rules for one- and 
two-gluon emission vertices are known for the NRQCD 
and AsqTad actions and are used in the present calcu- 
lations as well, the HISQ vertices needed to be handled 
differently. Analytic expressions for HISQ vertices are 
too complicated to write down in closed form. Instead 
we build up one- and two-gluon emission vertices emerg- 
ing from the HISQ action from vertices of simpler oper- 
ators through repeated use of convolution rules [36 . For 
instance, since one- and two-gluon emission vertices are 
known for once fattened links from the AsqTad Feynman 
rules, one can use them to build up vertices of a prod- 
uct of three, five, seven such fat links and implement the 
second fattening. 

We use Mathematica to carry out all the Dirac al- 
gebra and also to take derivatives of NRQCD vertices 
with respect to external momenta. We have devel- 
oped Fortran "automatic differentiation" routines to take 
derivatives of HISQ vertices. The same bookkeeping used 
for repeated application of convolution rules allows us 
here to apply the chain rule of differentiation each time 
two expressions are multiplied and build up derivatives 
of the complicated HISQ vertices. 

In our second method the correct infrared singularities 
were isolated and in many cases handled with subtraction 
functions. Details of the subtraction functions are given 
in Appendix [B} 

We believe that these two methods are sufficiently in- 
dependent that, in conjunction with tests of gauge in- 
variance and correct infrared behaviour and the replica- 
tion of results in the literature, agreement between these 
methods provides a stringent check of our results. 

We now give our numerical results for the HISQ quark 
mass and wavefunction renormalisation. 



1. Massless Quarks 

For massless quarks we require only the wavefunction 
renormalisation. In this case amtree = etrcc = and so 



Equation ( 53 ) reduces to 



d(apo) 



zsin(apo)4^^ +SW 



(55) 



The wavefunction renormalisation is infrared divergent 
and we decompose our results into an infrared finite con- 
tribution, C„, and an infrared divergent contribution. 



TABLE II. One-loop wavefunction renormalisation for 
massless HISQ fermions. The gauge parameter is ^. All 
uncertainties are statistical errors arising from the numerical 
integration of the relevant diagrams. 





y^rbow 




c. 


1 


-0.8183(1) 


0.4243(3) 


-0.3940(3) 





-0.0198(1) 


0.1343(3) 


0.1145(3) 



C^^ . Thus we write 



Zq = 1 

= 1 



a,zW+0(«2) 

as (Cf + +0(a2). 



The infrared divergence is given by 



IR 



1 

3^ 



[l + (e-l)]log (a^A^) 



(56) 



(57) 



where A is the gluon mass, introduced to regulate the 
infrared behaviour, and ^ is the gauge parameter. For 
massless quarks the infrared divergences in the lattice 
matching coefficients, arising from the wavefunction and 
vertex renormalisations, are ultimately cancelled by cor- 
responding divergences in continuum QCD. We confirm 
that any gluon mass dependence cancels between the lat- 
tice and continuum one loop coefficients. 

In contrast to the AsqTad and NRQCD actions, we do 
not need to use tadpole improvement for HISQ and the 
only contributions to the infrared finite piece, Cq, are the 
rainbow and tadpole diagrams. 



rbow 



C 



'tad 



(58) 



We tabulate our results for the wavefunction renormal- 
isation in Table HU 



2. Massive Quarks 

We require both the mass and wavefunction renormal- 
isation for massive HISQ fermions. In general, both of 
these quantities are functions of ei. For consistency with 
the HISQ action used in numerical simulations, however, 
we ignore ei in Equations (46 1 and (53). In Reference 



jl3l it was found that the nonperturbatively determined 
values for e were always close to etrco- This justified ig- 
noring one-loop (or higher order) corrections to etroo in 
all subsequent numerical simulations with massive HISQ 
quarks. Perturbative matching that is going to be com- 
bined with numerical computations must be consistent 
with the latter. We set e = etree accordingly. 

Neglecting ei considerably simplifies the perturbative 
calculation of both ami and Zq\ For completeness we 
tabulate our results for etroe, amtrcc and ami in Table 
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TABLE III. Tree level and one loop tuning parameters for 
massive HISQ fermions. AH uncertainties are statistical 
errors arising from the numerical integration of the relevant 
diagrams. 



amo 


^trcc 


amtroc 


ami 


0.826 


-0.344960900736 


0.814526131431 


0.6580(1) 


0.818 


-0.340115648115 


0.807017346575 


0.6551(1) 


0.645 


-0.234829780198 


0.641330413102 


0.5871(1) 


0.6300 


-0.225853340666 


0.626715862647 


0.5811(1) 


0.627 


-0.224064962178 


0.623789107649 


0.5795(1) 


0.6235 


-0.221981631663 


0.620372982565 


0.5784(1) 


0.6207 


-0.220317446966 


0.617638873348 


0.5771(1) 


0.434 


-0.117189612523 


0.433453860575 


0.4855(1) 


0.4130 


-0.106941294689 


0.412571424109 


0.4734(1) 


0.4120 


-0.106461983347 


0.411576478677 


0.4728(1) 



III We present results for a range of quark masses cor- 
responding to the MILC ensembles used in [6], [21] and 

m- 

The one loop mass renormalisation is gauge invariant 
and infrared finite, whilst the wavefunction renormalisa- 
tion is gauge dependent and infrared divergent. We write 
the one loop wavefunction renormalisation in Equation 
(531 as 



(59) 



where 



Zq^^^ = —Zq^ sinh(amtr 



(1 + etroc) 



{2 (cosh( 



a'^trco)) + (sinh(amtreo))^ 



d(apo) 



sin(apo)So^'' 



(I) 



(60) 
(61) 



Recall that we have set ei = 0. The contribution from 
Zq contains the logarithmic infrared divergence. In line 
with our presentation of the massless case, we separate 
the infrared finite and divergent pieces of the one loop 
self energy-dependent contribution, which we denote Cq 
and C^Q respectively. Thus we have 



.(1) _ 



c, 



IR 



(62) 



where the infrared divergent contribution is given by 



= ^ [-2 + (e - 1)] log (a^A^) 



(63) 



We further decompose the infrared finite contribution 
into the self energy rainbow and tadpole diagram and 
TOi-dependent pieces: 



Cq = ZX 'ami 



Cq + Cq . (64) 
We give our results for the one loop wavefunction 



renormalisation in Feynman gauge in Table IV 



renormalisation for 
in Feynman gauge. 



TABLE IV. One-loop wavefunction 
massive HISQ fermions. All results 
Quoted uncertainties are statistical errors from the numerical 
integration of the relevant diagrams. The quantity Z'q^^\ 
defined in Equations (60 1 and 1 47 1, is an analytic function of 
only etroc and amtrco- These parameters, given in Table [Till 
are known to twelve significant figures; we therefore neglect 
the uncertainty in Zq^^^ here. 



amo 


y^rbow 




7(mi) 


Cq 


0.826 


-1.342(1) 


0.1952(1) 


0.427495 


-0.865(1) 


0.818 


-1.349(1) 


0.1989(1) 


0.415945 


-0.878(1) 


0.645 


-1.510(1) 


0.2888(1) 


0.210922 


-1.097(1) 


0.6300 


-1.511(1) 


0.2949(1) 


0.197029 


-1.102(1) 


0.627 


-1.530(1) 


0.2970(1) 


0.194322 


-1.120(1) 


0.6235 


-1.534(1) 


0.2981(1) 


0.191192 


-1.125(1) 


0.6207 


-1.537(1) 


0.2982(1) 


0.188712 


-1.130(1) 


0.434 


-1.785(1) 


0.3652(1) 


0.066076 


-1.388(1) 


0.4130 


-1.820(1) 


0.3715(1) 


0.057058 


-1.421(1) 


0.4120 


-1.821(1) 


0.3712(1) 


0.056650 


-1.423(1) 



E. NRQCD Parameters 

The one loop parameters of NRQCD have been exten- 
sively studied in the literature, for example in pTl[^[5S] . 
Indeed, a two loop calculation of the energy shift has re- 
cently been carried out with a mixed approach that com- 
bines automated lattice perturbation theory calculations 
of the fermionic contributions with results extracted from 
quenched weak coupling simulations for all other con- 
tributions [37] . Here we simply introduce the notation 
and summarise the necessary results at the heavy quark 
masses relevant for the simulations in fSI . We require the 
wavefunction renormalisation, Zh, and the mass renor- 
malisation, Zm'- 



Zh = 1 
Zm - 1 



a, {CYt + Ch) 



(65) 
(66) 



The infrared behaviour of NRQCD must match that of 
full QCD and is therefore just 



1 



-2 + (e - 1)] log (a^A^ 



(67) 



In this case the infrared finite contribution, Ch, is 
composed solely of the heavy quark rainbow diagram, 
because both the tadpole diagram and tadpole improve- 
ment contribution vanish [22 . The mass renormalisation, 
on the other hand, depends on both the rainbow and tad- 
pole diagrams and the tadpole improvement term, 



C 



M 



Crbow 
M 



Ctad 
M 



/-mo 



(68) 
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TABLE V. One-loop heavy quark parameters in NRQCD. 
All results use stability parameter n — 4. We implement 
tadpole improvement with the Landau link definition of uo- 
All results are in Feynman gauge. The quoted uncertainties 
are statistical errors from the numerical integration of the 
relevant diagrams. 



aMo 


Ch 


Cm 


3.297 


-0.235(1) 


0.167(1) 


3.263 


-0.241(1) 


0.176(1) 


3.25 


-0.244(1) 


0.178(1) 


2.688 


-0.362(1) 


0.262(1) 


2.66 


-0.366(1) 


0.264(1) 


2.650 


-0.371(1) 


0.267(1) 


2.62 


-0.374(1) 


0.272(1) 


1.91 


-0.617(1) 


0.434(1) 


1.89 


-0.627(1) 


0.448(1) 


1.832 


-0.657(1) 


0.466(1) 


1.826 


-0.660(1) 


0.468(1) 



where an analytic expression for is given in 
3 



- 1 



3(aMo) 



2n{aMo) 
f '^1 



C5 

3 



\{aMor 



2n(aA/o)2 



,(1) 



(69) 



At one loop we need not distinguish between the pole 
mass and the bare mass, so for convenience we express 
all results in terms of the bare mass. 

We tabulate our results for Ch and Cm in Table [Vl 
We present results with Ci = 1 and use the Landau link 
definition of the tadpole improvement factor ug, with 



,(1) 



0.7503(1). All results use stability parameter 



IV. THE MATCHING PROCEDURE 

In lattice QCD the axial-vector and vector current op- 
erators mix with higher order operators under renormali- 
sation. In this section we outline the perturbative match- 
ing procedure that relates the lattice and continuum cur- 
rents and the extraction of the one loop mixing matrix 
elements. 

Our strategy for the perturbative matching of heavy- 
light currents with massless relativistic quarks and 
NRQCD heavy quarks follows that developed in 
and outlined in [32]. We will briefly review the match- 
ing formalism and refer the reader to the earlier articles. 
A related matching calculation for massless HISQ quarks 
with NRQCD formulated in a moving frame (mNRQCD) 
was undertaken for the vector and tensor heavy-light cur- 
rents in [31] . 

For massive quarks similar matching calculations using 
the same lattice action for both quarks have been carried 



out for Wilson quarks in f35] and for various implemen- 
tations of NRQCD in t30n40l[42]. To our knowledge, no 
matching calculations with mixed actions and massive 
relativistic quarks have been reported in the literature. 

Moving from massless to massive relativistic quarks 
complicates the matching procedure. In the former case, 
quarks and antiquarks at zero spatial momentum are in- 
distinguishable and consequently scattering and annihila- 
tion processes give identical results. In the massive case, 
however, we must distinguish between quarks and anti- 
quarks. For HISQ quarks at zero spatial momentum, this 
corresponds to choosing apo = ii'mtiee or apo = —iarritree 
respectively. We choose outgoing quarks or antiquarks — 
the "scattering" or "annihilation" channels respectively 
— to ensure we do not attempt to compute vanishing 
matrix elements. Thus we calculate the matrix elements 
of Vq and Ak in the scattering channel and of Aq and 
Vjt in the annihilation channel. This procedure is valid, 
even at non-zero lattice spacing, provided we match to 
the same channel in continuum QCD. 

Unfortunately, from the practical viewpoint of calcu- 
lating Feynman diagrams, using massive quarks compli- 
cates the numerical integration considerably. The chief 
difficulty lies in the annihilation channel, which contains 
a Coulomb singularity that must be handled with a sub- 
traction function. We discuss the subtraction functions 
employed in this work in more detail in Appendix |B] 
Furthermore, in the automated perturbation theory rou- 
tines, the pole in NRQCD propagator crosses the con- 
tour of integration and we can no longer carry out the 
usual Wick rotation back to Minkowski space. Instead, 
we must deform the integration contours and introduce 
a triple contour to ensure the stability of numerical inte- 
gration [501 IHT] . 

For both channels, the lattice matrix elements must be 
matched to their continuum QCD counterparts. Analytic 
expressions for the relevant QCD contributions already 
exist in the literature. References |;40 and [41 discuss the 
annihilation channel for the axial- vector current, whilst 
[32] present results for both vector and axial-vector cur- 
rents in the scattering channel at non-zero spatial mo- 
mentum. Results for components of both currents at 
zero spatial momentum in both channels are presented 
in [3S]. Whilst the authors of [30] are also concerned 
with calculating matching coefficients for the spacelike 
components of the vector current for lattice NRQCD, a 
procedure conceptually similar to that discussed in this 
work, they take a slightly different approach, calculating 
the continuum integrals numerically 

We calculate the mixing matrix elements required to 
match the axial-vector and vector currents in the effec- 
tive NRQCD theory to full QCD for the following com- 
binations of currents, Lorentz indices and orders in the 
perturbative and I/Mq expansions; 

1. massless relativistic quarks: 

(a) Vq through 

©(as, Aqcd as/(aMo), asAqcD/Mo); 
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(b) Vfe (fc = 1,2,3) through 

0(a.,AQCD/Mo,as/(aMo)); 

2. massive relativistic quarks: 

{a) Vf, {i^ = l,...,4:) through 
0{as, Aqcd/ Mo, as /{aMo)); 

(b) Af, through C'(a„AQCD/Mo,a,/(aMo)). 

The results for both axial-vector and vector currents are 
identical for massless relativistic quarks. To simplify our 
presentation we therefore only give results for the vector 
current for massless HISQ quarks. 

We discuss each of these cases in turn. 



A. Massless Quarks 

1. Temporal vector current 

We require three currents to match the temporal 
component of the vector current on the lattice to full 
QCD through ©(a,, Aqcd/Mq, a./(aMo), a.AqcD/Mo)- 
These are 



jj^\x)=q{x)r^Q{x), 



^,(.)r,,.^0(x). 



(70) 
(71) 



(^) = - ^Pio) «(^)^ • ^7or^Q(2:). (72) 

Here the Q fields are four component Dirac spinors with 
the upper two components given by the two component 
NRQCD field and lower components equal to zero. The 
r^t operator represents the vector current operator, so 
that here we have = TV- The difference operator V 
is defined in Equation ( |11[ ), with the arrow indicating 
whether the operator acts to the left or right. The Eu- 
clidean gamma matrices obey 

{7m'7i.} = 2(5^1/, 7^=7^- (73) 

The matrix element of the timelike vector current in 
full QCD is related to the matrix elements of the currents 
in the effective theory via 

(^o> = (i + (^r > + (i + asp^) (4''^-'') 

+ a.pf«)(jf-"'^)-l-0(a^,A|cD/Mo^a2a.). 

(74) 

Here we have expressed the lattice currents in terms of 
the subtracted currents. 



(75) 



for i = 1, 2. The subtracted currents are more physical 
and have improved power law behaviour j44j . 



The matching coefficients are given by 



P6 



(Vn) _ r(Vo) 



1 



(Vb) 
00 ' 



(76) 



iVo) 



Pi 



iVo) ^ q{Vo) 



-Q 



(Vo) 



/■(Vo) 
Sl2 J 



(77) 
(78) 



where the Bi arise from the matrix elements in full QCD 
and are given by [55 1 [55] 



B, 



(Vo) 



- (\n{aMo) 



1 



4 



B. 



(Vo) 



(79) 
(80) 
(81) 



The renormalisation parameters Cq, Ch and Cj\/ are the 
one loop self energy corrections discussed in the previous 
sections. For convenience we have written the pole mass, 
which is common to both lattice and continuum theories, 
in terms of the bare quark mass. We must therefore 
include the one loop mass renormalisation that relates 
these two masses in the tree level contribution from X 



(1) 



The C,r°^ in Equations (|76|) to (|78|) are the one loop 



mixing matrix elements that arise from the mixing of the 
currents. These contributions are generated by the one 
loop diagrams in Figure [2] From left to right these are: 
the "vertex correction" diagram, the "light earlobe" dia- 
gram, the "heavy earlobe" diagram and the "vertex tad- 
pole" diagram. To extract the mixing matrix elements, 

(i) 

we insert one of the lattice currents, Jg , at the vertex 
and then project out the tree level expression {JQ^jtree- 
Thus, for example, Cio represents the projection of Jq^^ 



onto Jq"' and Cii the projection of Jq''' onto itself. 

Some of the mixing matrix elements are infrared diver- 
gent and, as for the wavefunction renormalisation con- 
tributions, we separate the infrared divergent and finite 
pieces. For example, we write 



7(1) 



where 



Coo — Coo + CoOi 



(82) 



(83) 



We confirm that all infrared divergences ultimately can- 
cel in the matching coefficients pi. Demonstrating that 
the matching coefficients are infrared finite is a nontrivial 
check of our results. 

The matrix element C02 includes a term that removes 
an 0{aas) discretisation error from 4°^ [221138]. Thus 
the matching procedure ensures that ©(a^/Mo) and 
0{aas) corrections are made at the same time. 
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FIG. 2. Contributions to the one loop mixing matrix ele- 
ments to match the vector and axial-vector currents in lat- 
tice NRQCD to full QCD. Clockwise from top left to lower 
left are the "vertex correction" diagram, the "heavy earlobe" 
diagram, "vertex tadpole" diagram and the "light earlobe" 
diagram. The double lines indicate heavy quarks, the single 
lines represent light quarks and the wavy lines are gluons. 
Operator insertions are denoted by the solid circles. 



Finally we note that there is a second dimension four 

(2) 

current operator that is equivalent to Jq via the equa- 



tions of motion 



J, 



•(2) 



1 ^ 



(84) 



where the arrow indicates that the derivative acts to the 
left. The effects of this current operator must be included 
in the determination of ^^2- 



2. Spatial vector current 

In this case we require only the first two of the three 
lattice currents given above, those of Equations ( 70 1 and 



(|71|). The matrix element of Vk in full QCD is related to 

(85) 



the effective NRQCD current via 



where 



and 



t(1), sub 



(14 



(14) 



11' 



(14 



/■(14) 
'^00 ' 



(86) 



(87) 



ln(aMo)-^2 



The only contribution to Cojl!"'' ^^"^ CicT*'' vertex 
correction diagram in Figure with the current or 
Jn^"* inserted at the vertex. 



B. Massive Quarks 

The matching calculation for massive HISQ quarks 
proceeds in a similar manner to the massless case just 



discussed. Here, however, one must rescale the lattice 
currents, J/i , by the tree-level massive HISQ wavefunc- 

tion renormalization ( Zq \ .In the following we as- 
sume that the currents have been rescaled. 



1. Vector current 



We again require only two of the three lattice currents: 
J/l°^ and . We write the matrix element of the vec- 

tor current in full QCD in terms of the matrix elements 



of Jq"^ and J^i"^^'^"'' as 



where, in this case. 



t(1), sub 



T(l),SUb _ t(1) 



(89) 



(90) 



We denote the matching coefficient for massive HISQ 
quarks by 770, to clearly distinguish the massless and mas- 
sive cases. The matching coefficient is given by 



(Vm) 



1 



with 



D, 



,(Vo) 


1 


f amtroc 


+ aMa 


In 


/ amtrcc 





TT 


\ amtrcc 


- aMo 


[ aMo 


(Vfc) 






- aMo 




^ amtrcc 







V O.ITI'trcc 


+ aMo 


V aMo 



(91) 

(92) 
(93) 



The Tij are the mixing matrix elements for massive rela- 
tivistic quarks. 

The leading order mixing matrix elements for the tem- 
poral vector current are logarithmically infrared diver- 
gent. Hence we write 



where 



:(Vo) 



TR 



.{Vo) 



IR 
'00 J 



'00 = ^log(a^A"). 



(94) 



(95) 



In contrast to the massless case the infrared divergences 
in the vertex and wavefunction renormalisations cancel 
separately in both the lattice and continuum matrix ele- 
ments. Confirming that the sum of the lattice results is 
infrared finite serves as an important crosscheck of our 
calculation. 

The evaluation of the mixing matrix elements for the 
spatial vector current is more complicated than for the 
temporal component. In this case the mixing matrix el- 
ement Tq^''^ contains not only a logarithmic divergence 
but a linear divergence as well: 

^IR 



'00 ~ '00 



'00 



1 



Stt 



'00 ' 
amtrccaMo 



1 



+ 21og(a2A2) ) . (96) 



amtree + oMq aX 
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The logarithmic divergence is cancelled by the wavefunc- 
tion renormalization, leaving both lattice and continuum 
contributions with a linear divergence. This, in turn, 
cancels when we match the lattice and continuum results 
so that the matching coefRcient is infrared finite. 

2. Axial-vector current 



TABLE VI. One-loop matching for the temporal vector 
current with NRQCD heavy quarks and massless HISQ light 
quarks. All results use stability parameter n = 4 in the 
NRQCD action. We implement tadpole improvement with 
the Landau link definition of uq. For the p^'^"' the quoted 
uncertainties are the errors from each contribution added in 
quadrature, whilst for Cio^"' uncertainty is purely the 
statistical error from numerical integration. 



The matching relation for the axial-vector current is 
given at leading order by 



Here we have 



where (35 



00 ' 



amtrcc + ciMq ( arrit 



flTOtrcc 



- aMo 
aMo 



TT V aTOtroc + oAfo 



aMt 



(97) 
(98) 

(99) 
(100) 



Here the Aq current develops a linear IR divergence, 
which is the same as that given in Equation ( 96 1 for 
Vfe. This divergence again cancels between lattice and 
continuum results. 

In the following section we present our results for the 
matching coefficients pi of Equations ( 74 1 and ( 85 1 and 



770 of (891 and (97 1, together with the mixing matrix 
elements Cio and rio needed to fix J^j^^ 



aMo 








a(Vo) 

'>10 


3.297 


-0.072(2) 


0.048(2) 


-1.108(4) 


-0.0958(1) 


3.263 


-0.075(2) 


0.046(2) 


-1.083(4) 


-0.0966(1) 


3.25 


-0.075(1) 


0.046(2) 


-1.074(4) 


-0.0970(1) 


2.688 


-0.109(2) 


0.013(2) 


-0.712(4) 


-0.1144(1) 


2.66 


-0.110(2) 


0.013(2) 


-0.698(4) 


-0.1156(1) 


2.650 


-0.112(2) 


0.013(2) 


-0.696(4) 


-0.1157(1) 


2.62 


-0.116(2) 


0.008(2) 


-0.690(4) 


-0.1171(1) 


1.91 


-0.161(2) 


-0.038(3) 


-0.325(4) 


-0.1539(1) 


1.89 


-0.162(2) 


-0.038(3) 


-0.318(4) 


-0.1553(1) 


1.832 


-0.162(2) 


-0.042(3) 


-0.314(4) 


-0.1593(2) 


1.826 


-0.163(3) 


-0.043(3) 


-0.311(4) 


-0.1595(2) 




FIG. 3. Tadpole correction contribution to the one loop 
mixing matrix element pi. The double lines indicate heavy 
quarks, the single line the light quark and the cross represents 
the tadpole insertion. 



V. MATCHING PROCEDURE RESULTS 

As we discussed in the previous results section for the 

we im- 



quark renormalisation parameters. Section HID 



plement two independent calculation procedures to cross 

check our results. We have calculated all the relevant 

V r 

mixing matrix elements, C,^^^ and T^f^ , required to match 

the lattice currents with continuum QCD. For clarity of 

presentation, however, we only give our results for the 

V r 

final matching coefficients, ^ and rj^ . We also in- 

V r 

elude the mixing matrix elements, Cio and t-^^\ for com- 
pleteness, because these are needed to construct the sub- 
tracted lattice currents Jf;""'^. 



Cii^"^ and is given by 



^11 



,(1) 



(101) 



We use the Landau link definition of the tadpole correc- 
tion coefficient, Ug^'' = 0.7503(1), when calculating pi. 



In Table VII we give our results for the matching coef- 
ficients for the spatial components of the heavy-light vec- 
tor current with massless HISQ light quarks and NRQCD 
heavy quarks. 



B. Massive Quarks 



A. Massless Quarks 



We tabulate our results for the matching coefficients 
, (Vb) ^;ff„,.„v,f ^occ^c T„Kl^ YJ 



Pi " at four different heavy quark masses in Table 
Only the matching coefficient pi has a tadpole correction 
coefficent, arising from the tadpole correction insertion 
illustrated in Figure [3j This correction contributes to 



In Table VIII we tabulate our results for the matching 
coefficients for Vq with massive HISQ light quarks and 
NRQCD heavy quarks. 

We give our results for the matching coefficients for Vk 
in Table |IX| And finally, in Tables |X] and |XI[ we present 
our results for the matching coefficients for ^0 and 
respectively. 
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TABLE VII. One-loop matching coefficients for the spatial 
vector current with massless HISQ light quarks. See the 
caption accompanying Table fVl] for more details. 



aMo 






SIO 


3.297 

3.263 

3.25 

2.688 

2.66 

2.650 

2.62 

1.91 

1.89 

1.832 

1.826 




-0.046(2) 
-0.045(2) 
-0.045(2) 
-0.034(2) 
-0.034(2) 
-0.034(2) 
-0.033(2) 
0.007(2) 
0.009(2) 
0.020(2) 
0.020(2) 


0.0319(1) 
0.0322(1) 
0.0323(1) 
0.0382(1) 
0.0385(1) 
0.0386(1) 
0.0391(1) 
0.0513(1) 
0.0518(1) 
0.0532(1) 
0.0532(1) 


TABLE VIII. One-loop matching coefficients for the temporal 
vector current with massive HISQ quarks. All results use 
stability parameter n = 4 in the NRQCD action. For ryg^"' 
the quoted uncertainties are the errors from each contribu- 
tion added in quadrature, whilst for t^^°^ the uncertainty is 
purely the statistical error from numerical integration. 


aMo 


amo 




^10 


3.297 
3.263 
2.688 
2.660 
2.650 
2.650 
2.620 
1.910 
1.832 
1.826 


0.8260 
0.8180 
0.6300 
0.6450 
0.6235 
0.6207 
0.6270 
0.4340 
0.4130 
0.4120 


-0.151(3) 
-0.148(3) 
-0.121(3) 
-0.117(3) 
-0.113(3) 
-0.112(3) 
-0.116(3) 
-0.102(3) 
-0.098(3) 
-0.098(3) 


-0.0488(1) 
-0.0494(1) 
-0.0647(1) 
-0.0648(1) 
-0.0658(1) 
-0.0659(1) 
-0.0663(1) 
-0.0990(1) 
-0.1043(1) 
-0.1046(1) 


TABLE IX. One-loop matching for spatial vector current 
with massive HISQ quarks. See the caption accompanying 
Table IVIIII for more details. 


aMo 


amo 






3.297 
3.263 
2.688 
2.660 
2.650 
2.650 
2.620 
1.910 
1.832 
1.826 


0.8260 
0.8180 
0.6300 
0.6450 
0.6235 
0.6207 
0.6270 
0.4340 
0.4130 
0.4120 


-0.124(5) 
-0.118(5) 
-0.025(5) 
-0.024(5) 
-0.015(5) 
-0.014(5) 
-0.019(5) 
0.049(5) 
0.059(5) 
0.060(5) 


0.0420(1) 
0.0423(1) 
0.0484(1) 
0.0488(1) 
0.0489(1) 
0.0489(1) 
0.0493(1) 
0.0618(1) 
0.0636(1) 
0.0638(1) 



TABLE X. One-loop matching for the temporal axial-vector 
current with massive HISQ quarks. See the caption accom- 
panying Table [VHl] for more details. 



aMo 


amo 






3.297 
3.263 
2.688 
2.660 
2.650 
2.650 
2.620 
1.910 
1.832 
1.826 


0.8260 
0.8180 
0.6300 
0.6450 
0.6235 
0.6207 
0.6270 
0.4340 
0.4130 
0.4120 


-0.237(5) 
-0.232(5) 
-0.188(5) 
-0.192(5) 
-0.183(5) 
-0.182(5) 
-0.189(5) 
-0.219(5) 
-0.222(5) 
-0.221(5) 


-0.1260(1) 
-0.1269(1) 
-0.1452(1) 
-0.1464(1) 
-0.1468(1) 
-0.1467(1) 
-0.1480(1) 
-0.1853(1) 
-0.1908(1) 
-0.1914(1) 


TABLE XI. One-loop matching for the spatial axial-vector 
current with massive HISQ quarks. See the caption accom- 
panying Table |VIII| for more details. 


aMo 


amo 


,(-^) 


^10 


3.297 
3.263 
2.688 
2.660 
2.650 
2.650 
2.620 
1.910 
1.832 
1.826 


0.8260 
0.8180 
0.6300 
0.6450 
0.6235 
0.6207 
0.6270 
0.4340 
0.4130 
0.4120 


-0.260(3) 
-0.260(3) 
-0.194(3) 
-0.191(3) 
-0.183(3) 
-0.182(3) 
-0.185(3) 
-0.091(3) 
-0.076(3) 
-0.076(3) 


0.0163(1) 
0.0165(1) 
0.0216(1) 
0.0216(1) 
0.0219(1) 
0.0320(1) 
0.0221(1) 
0.0330(1) 
0.0348(1) 
0.0349(1) 




VI. 


SUMMARY 





We have calculated the one loop matching coefhcients 
required to match the axial- vector and vector currents on 
the lattice to full QCD. We used the HISQ action, with 
both massless and massive quarks, for the light quarks 
and NRQCD for the heavy quarks. As part of the match- 
ing procedure we have presented one loop mass and wave- 
function renormalisations for both HISQ and NRQCD 
quarks. We find that the perturbative coefficients are 
well behaved and none are unduly large. 

The matching coefficients for HISQ-NRQCD currents 
with massless HISQ quarks are important ingredients in 
the determination of heavy-light mesonic decay parame- 
ters from lattice QCD studies [6^. Recent studies of the 
Bs meson using the relativistic HISQ action for both b 
and s quarks have been carried out . Such an approach 
has the advantage that perturbative matching, which is 
generally the dominant source of error in the extraction 
of decay constants, is not required. Currently, however. 
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simulations at the physical b quark mass are prohibitively 
expensive and an extrapolation up to the b quark mass 
is still needed. Furthermore, simulations of the B me- 
son are not presently feasible, as the use of light valence 
quarks and close-to-physical b quark masses requires both 
large lattices and fine lattice spacings. In light of these 
considerations, the use of an effective theory for heavy- 
light systems remains the most efficient method for pre- 
cise predictions of Jb^ / Jb and /s . Such calculations re- 
quire the perturbative matching calculation reported in 
this article. 

The matching calculations reported in this work are 
also crucial for the HPQCD collaboration's nonpertur- 
bative studies of the semileptonic decays of B and Bg 
mesons with NRQCD and HISQ quarks. On the one 
hand, matching coefficients with massless HISQ quarks 
are required for the determination of the B — > iriiy, 
B — )■ K£~^£~ and Bg K£i/ decay parameters |7j. On 
the other hand, our results for the matching coefficients 
with massive HISQ quarks will be needed in future cal- 
culations of the B — > Div and Bg — ?> Dsiiy decay param- 
eters. 



sinh(£;) - l±i^££^ (sinh(£;))' 



which, for notational convenience, we write as 
Using the relations 



(A2) 



(A3) 



Yip, - 0) - 0, 



dE 
dpz 



= 0, 



dY{p,) 



dpz 



= 1, 



and differentiating twice using 



d _ d dE d 

dpz dp, dpz dE ' 



(A4) 
(A5) 



we find 
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Appendix A: The HISQ tuning parameter 

In this appendix, we derive expressions for the tree 
level and one loop tuning parameters, etrcc and ei. 
Throughout this appendix we neglect factors of the lat- 
tice spacing, a, for convenience. 

For an onshell particle with momentum given by p^ — 
{iE,0,0,pz) one defines the kinetic mass as 



d^E 
dpi 



(Al) 



Pz=0 



At nonzero momentum the tree level pole condition be- 



-X- 



dX\ fd'^E 



dE J \ dpi 



dY 

dpx 



= 0, (A6) 



and thus the tree level kinetic mass is 



^dX 



m(°) -X- 

= cosh(mtreG) sinh(mtroo)(l - ©)(! - 39). (A7) 



Here we have defined (1 -I- etrcc) (sinh(mtreo))^/6. 

Requiring r7i[.j^ = Titrcc imposes a condition on thi 
tree level tuning parameter that leads to Equation ( 45 ) 



At one loop, the procedure is much the same. The one 
loop pole condition is 



(mo - a^S*"))^ = X{E) - sinh(£;)a^E^' 



where 



X{E) = siiih{E) 



Y{p^_) - sm{p^_)asl^ 



1 - ^ {sinh{E)f 



(A8) 



(A9) 



Differentiating twice using Equation (A5) leads, after some algebra, to 



:dX 



TOkin = ~ '^A ^0 [^Q J Sinh(mtrcc) + 2TOtrocS^^^ ~ ^°~dE 



sinh(i;)4^) - sw" 



TOomtrco 



dpi 



sinh(mtrcc)s[,^^ - S^'^ 



(AlO) 
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where we have only kept terms up to 0{as)- 
For convenience, we write this as 



~dX 



(All) 



Using the expansions of Equations (42), (43) and (46 1, we can evaluate the product, at one loop to obtain 



~dx 



where 



1 3 / 

ui = - cosh(mtioc)(sinh(rntreo)) (^(1 + etrcc)(sinh(mtroc))^ 

U2 = (sinh(mtroc))^(l - ^-^^^|^(sinh(mtroo))^) 



1 + etr 



2(cosh(TOtrec))^ + (sinh(mtrec))' 



+ (cosh(mtrco))^(l - 



1 + etrc 



-(sinh(TOtroc))' 



(A12) 



(A13) 



(AM) 



We can therefore write Equations (All) and (A12) as 



with i^*-^) and z/^*^' given by 



and 



' = Ui + U2'm\ , 



mf^ = sinh(mtrcc)4''^-Sm, 



1 



(A15) 

(A16) 
(A17) 

(A18) 
(A19) 



We obtain an expression for the one loop tuning pa- 
imeter by e( 
The result is 



rameter by equating the one loop masses: m^|^ = m^^\ 



(A20) 



where v'^'^^ and i/^'' are given in Equations (A16l and 



(A17) respectively. 



Appendix B: Subtraction functions for numerical 
integration 

At intermediate stages of the lattice-to-continuum 
matching procedure one encounters infrared (IR) di- 
vergent integrals and care is required to ensure that 
Vegas can handle them accurately. For diagrams involv- 
ing massless HISQ fermions, it is usually sufficient to 
introduce a nonzero gluon mass A, fit results to appro- 
priate functions of this mass and then extract the IR 



finite parts. For massive HISQ fermions, on the other 
hand, it is often necessary to include specific subtrac- 
tion terms into the integrand in order to stabilize the 
Vegas integrations. In this appendix we list examples of 
such subtraction terms. Given an IR divergent integral, 



(2^ •^iat(fc,A), 



(Bl) 



where Cp = 4/3 is a color factor (the quadratic Casimir 
operator) required to correctly normalize the infrared di- 
vergences. We employ subtraction terms in the following 
way: 



I=Cf J ^^^^4 {J"lat(fc,A) - J"sub(fc,"le//,A, A)} 



(27r) 

-f F(771eff,A,A), 



(B2) 



where 

f d'^k 

F(rnoff, A, A) = Cf J j^-^ J^suh{k,meff,A, X). (B3) 

Here A is a cutoff imposed on J^suh such that J^suh = for 
k^ > A^, and m^ff is defined below. The full expression 
for Z in (B2) does not, of course, depend on A. We have 
done the calculations with two different values for A, e.g. 
aA = 2 and 3, to check this. 

The choice for J-guh is far from unique. One wants a 
function of fc^ with the same IR behaviour as the original 
integrand J^iat, that is, however, simple enough that the 
integral in the "addback" function, F(r7ioff, A, A), can be 
evaluated with relative ease. One natural choice is the 
integrand of the corresponding continuum theory Feyn- 
man diagram J-"cont- This is what has often been done 
in the literature. For massive fermions there remains the 
question of what fermion mass to use in J^3ont■ It was 
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suggested in [5S] to pick a mass, denoted by nies, such 
that J-"sub(^, JTieff, A, A) mimics as closely as possible the 
correct fc^ — >■ limit in the denominator of the lattice 
fermion propagator. For instance in the continuum the- 
ory one would have (we work in Euclidean space), for 
onshell quarks with external momentum {im, 0), a 

fermion propagator with denominator given by 

denom ^ {p - kY + m'^ ^ - 2iniko. (B4) 

Taking a hint from ( |B4[ ) we pick nics by first setting the 
external momentum to — (imticoiO), expanding the 
denominator of the free HISQ propagator around kg = 
and then looking for the coefficient of (— 2zfco). One finds, 

rrioff = cosh(TOtroG) sinh(mtroc)(l - 6)(1 - 36), (B5) 

here Q is defined afte 
given in ( A7 ) , so that 



a result that may not come as a surprise. We note, how- 



where Q is de fined after (|A7). We recognize this as m^.;^. 



(0) _ 

"loff = ml^^ = mtree, 



ever, that the last equality in ( B6 ) holds only because we 



have tuned etrco to ensure 
was not the case for uses of nicff in the past 
volving massive Clover fermions. 



,(0) 



This 
51 in- 



Following the guidlines described above, the subtrac- 
tion term for the rainbow diagram correction to the mas- 
sive HISQ wave function renormalization Zq becomes 



^Z^ ^ ,(.2 _ ,2. r 4(fcg+bV4)((fc^)^-b^fc§) 
^ ^ 1 (A:2 + A2)((fc2)2 +52^2)2 

' fc2(fc2 + A2)((fc2)2+;,2fc2)J' V ) 



(B6) with b = 2mtrcc- This leads to an addback function 



2 

3^ 



log 

-1) 



A2 
A2 

log 



log 



where Wq = \/\P- + A2 and VFi ee [A -I- Wq\. 



Wi 



log 



A2 



Wi 



3A2) 
-(262 



A2 



AWq 



[h" - 6A2) 



.A2) 



64 

2^ AVKq 

64 



(36^ + 2A2 



(B8) 



Similarly, for the one-loop vertex correction for a scat- 
tering diagram involving Vq one has in Feynman gauge. 



-^sub ~ 



g(A2-fc2){fc2(fc2+52) _ _^(fc2)2| 

[(fc2)2 + 62A:2] [fc2 + (p/(2M))2] [F + A2] " 



(B9) 

And for the annihilation diagram involving Vk one has. 



-^sub " 



g(A2-fc2){fc2(fc2^^2) ^ _^(fc2)2| 
[(fc2)2 + h^kl] [fc2 + (P/(2M))2] [/c2 + A2] " 



(BIO) 

The only difference between (B9) and (B10[) is the rel- 



ative sign between the two terms in the numerator, i.e. 
the sign of the term linear in mtree- This is as it should 
be, since for annihilation one has an incoming anti-HISQ 
quark with the on-shell condition p^ (— imtrccO) re- 
placing the outgoing HISQ quark of the scattering pro- 
cess. The two terms in the numerator each lead to linear 
IR divergent results which cancel in the case of Vg leav- 
ing just a logarithmic IR divergence. For Vk one ends up 
with an expression with both linear and logarithmic IR 
divergent terms as is required. We have not attempted 
to integrate or J\^^ in closed form to obtain ana- 
lytic expressions for the addback functions and F^*" . 
Instead we reduced the integrals to ID integrals in the ra- 
dial variable < /c < A and used Vegas again to evaluate 
them. 
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